©AGU.PUBLICATIONS 


Journal of Geophysical Research: Planets 



RESEARCH ARTICLE 

10.1 002/201 3JE004502 


Key Points: 

• The zonal distribution of dust forcing is 
key to understanding wave response 

• The zonal wavenumber 1 pattern of 
dust forcing is the key mechanism 

• High-latitude wave response is 
coupled to that of the tropics 


Thermal tides during the 2001 Martian 
global-scale dust storm 

Scott D. Guzewich 1 , R. John Wilson 2 , Timothy H. McConnochie 3 , Anthony D. Toigo 4 , 

Donald J. Banfield 5 , and Michael D. Smith 1 

'NASA Goddard Spaceflight Center, Greenbelt, Maryland, USA, 2 Geophysical Fluid Dynamics Laboratory NOAA, Princeton, 
New Jersey, USA, department of Astronomy, University of Maryland, College Park, Maryland, USA, 4 Johns Hopkins 
University Applied Physics Laboratory, Laurel, Maryland, USA, 5 Center for Radiophysics and Space Research, Cornell 
University, Ithaca, New York, USA 


Correspondence to: 

S. D. Guzewich, 
scott.d.guzewich@nasa.gov 


Citation: 

Guzewich, S. D., R. J. Wilson, T. H. 
McConnochie, A. D. Toigo, D. J. Banfield, 
and M. D. Smith (2014), Thermal tides 
during the 2001 Martian global-scale 
dust storm, J. Geophys. Res. Planets , 7 79, 
506-519, doi:1 0.1 002/201 3JE004502. 

Received 14 AUG 2013 
Accepted 7 FEB 2014 
Accepted article online 12 FEB 2014 
Published online 6 MAR 2014 


Abstract The 2001 (Mars Year 25) global dust storm radically altered the dynamics of the Martian atmosphere. 
Using observations from the Thermal Emission Spectrometer onboard the Mars Global Surveyor spacecraft and 
MarsWRF general circulation model simulations, we examine the changes to thermal tides and planetary waves 
caused by the storm. We find that the extratropical diurnal migrating tide is dramatically enhanced during the 
storm, particularly in the southern hemisphere, reaching amplitudes of more than 20 K. The tropical diurnal 
migrating tide is weakened to almost undetectable levels. The diurnal Kelvin waves are also significantly 
weakened, particularly during the period of global expansion at L 5 = 200°-210°. In contrast, the westward 
propagating diurnal wavenumber 2 tide strengthens to 4-8 K at altitudes above 30 km. The wavenumber 1 
stationary wave reaches amplitudes of 10-12 K at 50°-70°N, far larger than is typically seen during this time of 
year. The phase of this stationary wave and the enhancement of the diurnal wavenumber 2 tide appear to be 
responses to the high-altitude westward propagating equatorial wavenumber 1 structure in dust mixing ratio 
observed during the storm in previous works. This work provides a global picture of dust storm wave dynamics 
that reveals the coupling between the tropics and high-latitude wave responses. We conclude that the zonal 
distribution of thermotidal forcing from atmospheric aerosol concentration is as important to understanding 
the atmospheric wave response as the total global mean aerosol optical depth. 


1. Introduction 

Global-scale dust storms are one of Mars' most enigmatic features. Despite the repeatability of many aspects 
of the Martian climate from year-to-year, global-scale dust storms do not occur annually and are greatly 
variable in duration and intensity [e.g., Liu et at., 2003; Smith, 2008]. Martian global-scale dust storms have 
been observed to occur during the second half of the Martian year in the seasons bracketing planetary 
perihelion at L s = 251° when greater insolation warms the atmosphere and enhances dust lifting [ Kahn etal., 
1 992]. Typical global average visible wavelength column dust optical depths of approximately 0.3-0.5 are 
observed early in the Martian year when the atmosphere is less dusty. In dustier portions of the year column 
dust optical depth is typically twice that value [Smith et al., 2001]. For the 2001 global-scale dust storm, the 
Thermal Emission Spectrometer (TES) onboard the Mars Global Surveyor (MGS) spacecraft retrieved infrared 
(9.3 pm) optical depths of 1.5-2 [Smith et at., 2002], which translates to approximately twice that in visible 
wavelengths [Wolff et al., 2006], 

In 2001, or Mars Year 25 (MY25) in the calendar of Clancy et al. [2000], a dust storm developed unusually early in 
the southern hemisphere (L s = 1 84°) just after the southern spring equinox near Hellas Basin [Smith et al., 2002; 
Strausberg etal., 2005; Cantor, 2007]. The storm expanded eastward as additional areas of dust lifting developed, 
and a thick dust cloud encircled the planet by L s = 1 97° [Cantor, 2007]. Dust reached heights of 60-80 km and 
took months to completely settle out of the atmosphere back to nominal levels [Clancy et al., 2010], 

The dust storm forced dramatic changes to atmospheric temperatures and structures. While the daytime 
surface temperatures cooled by 20-30 K due to decreased insolation reaching the surface, atmospheric 
temperatures rose by 20-40 K, depending on height [Smith et al., 2002; Smith, 2004]. Smith et al. [2002] also 
noticed a dramatic diurnal variability in TES temperatures (20-50 K) during the storm and attributed it to 
enhanced thermal tides. Dust forcing of thermal tides has previously been examined in observations and 
general circulation models (GCMs). Leovy and Zurek [1 979] noted significant enhancement of the migrating 
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Figure 1 . TES 2 P.M.dust mixing ratio (Ar/mb) at 60 km altitude for MY25 (a) 4 = 170°, (b) 4 = 180°, (c) 4 = 190°, (d) 4 = 200°, 
(e) 4 = 210°, and (f) 4 = 220°. 


diurnal and semidiurnal surface pressure tides from Viking Lander pressure sensor data during both global- 
scale dust storms of 1 977. Haberle et al. [1 982], Wilson and Hamilton [1 996], and Wilson and Richardson [2000] 
later modeled the 1977 global dust storms and examined the role of tides. Indeed, some studies have linked 
tides to the development and growth of global-scale dust storms [Leovy and Zurek, 1979; Leovy, 1981; Basil 
etal., 2006; Wilson et al., 2008a, 2008b] and large regional storms [Wang etal., 2003]. Building on those earlier 
studies, we find that the zonal distribution of the dust forcing is as important to the resulting wave generation 
as is the total dust optical depth. This zonal distribution of dust is plotted at high altitude (60 km) in Figure 1 
from before the storm (L s = 170°) through the initiation of storm decay (4 = 220°). The slowly westward 
propagating dust maximum in Figure Id-lf was first identified by Clancy etal. [2010], which is also described 
in Guzewich et al. [2013a] and appears to be a key player in the atmospheric wave response during the storm. 

In this paper, we examine the effect of the 2001 global-scale dust storm on atmospheric thermal tides using 
TES observations and the MarsWRF GCM. In section 2 we describe the TES data set and the MarsWRF GCM. In 
section 3 we analyze and discuss the data and present our results. Finally, section 4 concludes the work. 

2. Data and Model Description 

2.1. TES Limb-Sounding Data Description 

TES retrieved temperatures through both nadir- and limb-sounding observing modes during its mission. This 
study uses primarily limb-sounding data, but nadir-derived data are incorporated into the retrieved 
temperature profiles. Indeed, below 80 Pa (approximately 22 km altitude), the data are 100% derived from 
nadir-viewing observation. The limb-viewing observations are then phased into the nadir observations 
above 80 Pa and above 10 Pa, all data are derived from limb observations. The primary difference to the TES 
temperature retrievals described by Smith et al. [2001 ] is that errors have been significantly reduced through 
accurate instrument background measurements and noise covariance matrices [McConnochie and Smith, 
2011] (for additional details, see Guzewich et al. [2013a, Appendix 1]). The temperature retrieval algorithm is 
performed first and the aerosol opacities retrieved second. At very high opacities, as was seen during the 
height of the 2001 global-scale dust storm, this may create larger temperature errors (the errors could be as 
large as 5-20 K) in the nadir-viewing geometry in the lowest 1-1.5 scale heights of the atmosphere (with 
errors increasing toward the surface). We have performed sensitivity tests of this scenario and find that above 
1.5 scale heights, the temperature errors are comparable to the low-opacity errors: 1-3 K, which is 2% or less 
of typical temperatures. Altitudes below 1.5 scale heights above the surface are not plotted in the figures. The 
limb-derived temperatures may also be biased near 4 scale heights (the lowest altitude where the 
temperature profile is almost completely derived from limb-viewing observations) through dust absorption 
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m = 0 


Table 1 . Waves Discussed and scattering from the surface. However, 

Field Type s a Comments at 5 scale heights and above, the dust 

opacity is too low to bias the temperature 
retrieval. The smooth continuity in 
temperatures and wave amplitudes 
between the two altitude regions of 
highest reliability (1.5-4 and above 5 scale 
heights) and corroboration from the 
general circulation model gives us high 
confidence in our results and interpretations. 

Data were binned by 10° latitude, 30° 
longitude, and 1 0° of solar longitude 
resulting in approximately 50 temperature 
profiles averaged together to produce one 
bin value. Because MGS flew a Sun- 
synchronous orbit, observations were 
constrained to be at approximately 2:00 A.M. 
and 2:00 P.M. local solar time. Our binning 
divided temperatures into either "morning" 
or "afternoon" (r 2am and r 2pm ). At very high 
latitudes, the local time of the observations deviates somewhat from the nominal 2:00 A.M. and 2:00 P.M. times. 
This may complicate our analysis at those latitudes, but at middle and low latitudes where our analysis is focused, 
the typical deviations from the nominal local times are less than 1 h and have negligible effect on this analysis. 
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In our analysis, we assume an atmospheric wave has functional dependence on longitude (1), pressure (p), 
latitude (4»), and time (f). In a Sun-synchronous frame of reference, local time is fixed. We neglect the pressure 
dependence to produce the following wave equation: 

T(X, fLsr) = lAya cos[(s - a)X + af LST + 4] (1) 


where A is the amplitude of a wave with wavenumbers and frequency per sol a, f LST is local solar time, and S is 
an arbitrary phase. Positive (negative) a denotes westward (eastward) wave propagation. Analysis of the tides 
proceeded as many previous studies [Banfield et al., 2003; Lee et al., 2009; Guzewich et al., 201 2] by creating 
"average" and "difference" temperatures of the morning and afternoon observations. 

f~avg — (f"2pm “F T 2am )/2 (2) 

Tdif — (f 2pm — T 2am )/2 (3) 


Based on the offset of n radians between the 2:00 P.M. and 2:00 A.M. local solar times, equations (1 )— (3) can be 
combined and solved for T avg and T di f, the solutions of which differ by only a factor of (— 1) CT . Thus, only terms 
with even (odd) a appear in the T avg (T djf ) fields (see Guzewich et al. [201 2] for additional details). For example, 
stationary waves (o = 0) and semidiurnal tides (o=2) appear in the T avg fields, while diurnal tides (a= 1) 
appear in the T dif fields. We then apply a Fourier transform to the T-average and T-difference temperature 
fields. The Sun-synchronous viewpoint aliases tides (and stationary planetary waves) into satellite-relative 
wavenumber fields (denoted by m henceforth). A wave observed in wavenumber m field can be any 
combination of waves where m = \s~o\ [Forbes and Hagan, 2000; Wilson, 2000; Lee et al., 2009]. Table 1 
identifies the waves discussed in this work and which wavenumber m field they manifest. In part, we attribute 
temperature amplitudes in a given T avg or T di f field to a given wave (tide or stationary wave) by the difference 
in phase progression direction in height (eastward or westward) and the rate at which the phase changes 
with height (slowly for a long wavelength, quickly for a short). 

2.2. TES T 1S Data Description 

In an effort to expand and corroborate the limb-sounding data set described in section 2.1, we utilize TES 
15 pm brightness temperatures obtained through nadir-viewing observations, referred to henceforth as T 15 . 
These data provide a depth-weighted view of atmospheric temperature, approximately representative of the 
50 Pa (25 km) level in the atmosphere [Wilson and Richardson, 2000]. They have longitudinal and temporal 
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Figure 2. (a) Zonally averaged TES T 15 T-average temperatures and (b) T-difference temperatures during MY25. 


resolution of 3° and 0.5°, respectively, both of which are much higher than the limb-sounding data set. Analysis 
of the T 15 proceeded exactly as that of the limb-sounding data with calculating average and difference 
temperatures and then taking a Fourier transform of those. The depth-weighting will provide an underestimate 
of amplitude of waves with short vertical wavelengths (e.g., the westward propagating diurnal nonmigrating 
tides) while accurately representing those with longer vertical wavelengths (e.g., the diurnal Kelvin waves). 

2.3. MarsWRF Description 

MarsWRF is the Mars-specific implementation of PlanetWRF, a generalized planetary atmosphere general 
circulation model developed from the terrestrial Weather Research and Forecasting (WRF) model [Richardson 
et al., 2007; Toigo et al., 201 2]. Cuzewich et al. [201 3b] describes how the dust distribution of MarsWRF was 
modified to utilize a dust climatology derived from the TES limb-geometry observations described by 
Guzewich et al. [201 3a]. For this work, we utilize a 3 Mars year simulation that "replays" Mars Years 24-27 and 
includes the 2001 global-scale dust storm. This is simulation S5 in the nomenclature of Guzewich et al. 
[2013b], and we will refer to it as the "dust storm simulation." We use a "standard" MarsWRF simulation 
(described as simulation SI by Guzewich etal. [2013b] and the "standard resolution" simulation by Toigo etal. 
[2012]) to serve as a control to compare against. The control simulation uses the "MGS scenario" Mars Climate 
Database dust prescription described by Montmessin et al. [2004], which uses a modified Conrath-v vertical 
dust profile [Conrath, 1975; Montmessin et al., 2004] and an analytic function describing the latitudinally and 
seasonally varying dust opacity typical of a year without a global-scale dust storm. Additionally, a more 
generic "dust storm" simulation is conducted where the vertical dust prescription of the "MGS Scenario" in 
the Mars Climate Database is used, but the visible band optical depth at a reference pressure of 610 Pa is 
globally set equal to 2. 

3. Analysis and Discussion 

3.1. Migrating Diurnal Tide 

The migrating diurnal tide is the westward propagating Sun-synchronous wavenumber 1 response to solar 
heating that occurs in planetary atmospheres [Chapman and Lindzen, 1 970]. As with Earth, the gravest Plough 
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Figure 3. Amplitude of the observed limb-sounding TES m = 0 T-difference field for (a) MY24 4 = 180°, (b) MY24 4 = 200°, 
(c) MY25 4=190°, (d) MY25 L s = 200°, (e) MY25 4 = 210°, and (f) MY25 4 = 240°. This represents the migrating diurnal tide. 


mode of the tide is equatorially trapped in the low latitudes on Mars (due to the apparent singularity near 30° 
where the Coriolis frequency equals the diurnal frequency). At higher latitudes, vertically trapped modes 
dominate. Previous modeling studies [e.g., Wilson and Richardson, 2000] have indicated that while the gravest 
Hough mode of the migrating diurnal tide is weakened by significant dust loading, the vertically trapped 
extratropical modes are strongly forced by atmospheric dust. 

In Figures 2a and 2b, we show the TES T 15 T-average and T-difference fields, respectively. The T-average 
temperature represents the zonal mean temperature and shows the rapid 20 K temperature increase in the 
lower atmosphere associated with the storm. The T-difference temperature demonstrates the intense 
thermotidal forcing associated with the storm from L s = 190° and 4 = 220° when optical depths were largest, 
causing the very large amplitudes in the middle latitudes of both hemispheres. Due to the nature of the T 15 
observations, the tropical migrating diurnal tide is not well represented (even prior to the storm), while the 
amplitudes of the vertically trapped extratropical modes are representative. 

However, as seen in the TES limb-sounding observations, the tropical migrating diurnal tide was weakened 
during the storm. In Figures 3a and 3b, we illustrate the m = 0 T-difference field of the TES limb-sounding 
observations at L s = 1 80° and 4 = 200°, respectively, for MY24. This provides a view of the "expected" nondust 
storm pattern to the migrating diurnal tide with a look of alternating minima and maxima in the vertical with 
amplitudes of approximately 8 K. A similar pattern is evident in MY25 at 4 = 180°, but the data coverage is 
significantly less than in MY24, and hence is not shown. This pattern is most apparent during the equinox 
seasons when the tropical migrating diurnal tide response is strongest. Similar structures have been 
observed with TES by Banfield et al. [2003] and Mars Climate Sounder by Lee et al. [2009] and Guzewich et al. 
[2012]. Immediately following the onset of the storm at L s = 190° in MY25 (Figure 3c), the extratropical tide 
becomes significantly stronger with amplitudes of 15 K in the southern hemisphere at 50°-60°S and 20-30 km 
and 10-12 K at a comparable location in the northern hemisphere. The equatorial structure is still present but 
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Figure 4. TES observed limb-sounding equatorial (10°S-10°N) temperatures anomalies (temperature minus zonal mean 
temperature) for MY25: (a) 2 A.M. L s = 200°, (b) 2 P.M. 4 = 200°, (c) 2 A.M. 4 = 210°, and (d) 2 P.M. 4 = 210°. 



has been weakened relative to 4 = 180°. As the storm became global near 4 = 200°, the observations indicate 
the vertically propagating tropical mode becomes completely suppressed (Figure 3d) and the extratropical 
vertically trapped modes are dominant through 4 = 290°. The amplitude of the extratropical tide peaks at 
4 = 200°-210° (Figures 3d and 3e) at 20-22 K in the southern hemisphere when globally averaged opacity was 
greatest, while the northern hemisphere (within latitudes of retrieved data) amplitude remains steady near 1 0- 
12 K. The observed southern hemisphere extratropical tide amplitude decreases soon after the storm reached 
maximum optical depth and dust lifting ceased (Figures 2b and 3f). Still, the extratropical tide amplitude 
remains anomalously high relative to nondust storm conditions, even for the dusty near-perihelion season, 
through northern fall and winter (compare Figure 3f with Guzewich et at. [2012, Figure 4c]). The vertically 
propagating tropical tide mode does not clearly reemerge until after 4 = 270°. 

3.2. Nonmigrating Tides 

Nonmigrating tides are those that are not Sun-synchronous, can propagate eastward or westward 
(depending on the particular tide), and are forced by a suite of atmospheric processes including interactions 
of the migrating tides and topography [Chapman and Lindzen, 1 970], tide-Rossby wave interactions 
[Moudden and Forbes, 2011a, 201 1 b], and others. 

Among the nonmigrating tides on Mars, the diurnal Kelvin wave 1 (DK1) has been the most thoroughly 
studied in the literature. Zurek and Leovy [1 981 ], Wilson and Hamilton [1 996], and Bridger and Murphy [1 998] 
examined the behavior of DKI's surface pressure signal in Viking Lander data during the 1977 global dust 
storms, and they found that the DK1 pressure tide was enhanced during the storm, leading to destructive 
interference between it and the westward propagating migrating diurnal tide. Martinez-Alvarado et at. [2009] 
and Montabone et at. [2008] (using a Mars data assimilation data set based on TES observations [Montabone 
et al., 2006]) speculated that a shift in phase of combined DK1 and migrating diurnal tide surface pressure 
signals, associated with dust heating near FHellas Basin from the onset of the 2001 storm, enhanced dust 
lifting in the Tharsis region and allowed the storm to strengthen. Still, the temperature amplitude of DK1 is 
expected to remain relatively invariant to zonally uniform aerosol loading due to a balancing between a near- 
resonant response in the cool dust-free atmosphere and stronger migrating tide forcing during enhanced 
aerosol loading [Zurek, 1988; Wilson and Hamilton, 1996; Guzewich et al., 2013b]. 

DK1 has been identified in many observations of the Martian atmosphere [e.g., Banfield et al., 2003; Wilson, 
2000, 2002; Hinson et al., 2001 ; Withers et al., 201 1 ; Guzewich et al., 2012] and exhibits amplitudes of 1 -4 K with 
a very long vertical wavelength resulting in no-phase progression with height (i.e., a barotropic structure). 
Due to the aliasing associated with a Sun-synchronous orbit, DK1 appears in the m = 2 T-difference field and 
will appear as a wavenumber 2 structure in longitude in the T 2am and T 2pm fields. This wavenumber 2 
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Figure 5. (a) MY25 4 = 210° amplitude of the observed limb-sounding TES m = 1 T-difference field, (b) diurnal westward 
wavenumber 2 tide amplitude from the MarsWRF global dust storm simulation, and (c) diurnal westward wavenumber 2 
tide amplitude from the MarsWRF control simulation, (d-f) The same fields for MY25 4 = 240°. 


structure is clearly apparent in the TES temperature fields in nearly all seasons with a sign change between 
2 A.M. and 2 P.M. (i.e., out of phase with each other, indicating a diurnal wave). During the 2001 global-scale 
dust storm, the dominance of DK1 over the nonmigrating tide temperature structure broke down. At L s = 1 90° 
in the TES observations, the equatorial (10°S-10°N) T 2am and T 2pm fields begin to lose the barotropic 
wavenumber 2 structure previously present and develop a wavenumber 1 structure with a strong westward 
phase tilt with height (not shown). The westward propagation in phase is indicated by the increasingly 
westward location of temperature maxima or minima with increasing altitude. Lewis and Barker [2005] and 
Wilson et al. [2008a, 2008b] found that the amplitude of DK1 in the surface pressure signal dropped during 
the 2001 global-scale dust storm (well after storm onset) in a Mars data assimilation data set described in 
Montabone et al. [2006]. Figure 4 shows the limb-sounding observed 2 A.M. and 2 P.M. temperatures (with the 
zonal mean temperature subtracted) for L s = 200°-210° (averaged between 10°S and 10°N) that most clearly 
shows this westward propagating wavenumber 1 field, particularly above 30 km altitude. We interpret this as 
a diurnal westward propagating wavenumber 2 nonmigrating tide (DW2), which appears as a wavenumber 1 
structure in a Sun-synchronous (fixed local time) reference and contributes to the m = 1 T-difference field. 

The gravest Plough mode of DW2 has a relatively short theoretical vertical wavelength of 3 1 .5 km and is 
equatorially confined to roughly 25° in each hemisphere (calculated using tidal theory as described in 
Chapman and Lindzen [1970]). Higher-order modes contribute to this tide in the extratropics. The m= 1 T- 
difference field shows the strong amplification of this tide with 4-8 K amplitudes at higher altitudes (above 
30 km) at E s = 200°-210° (Figure 5a). Also contributing to the m = 1 T-difference field is the zonally symmetric 
diurnal tide (see Table 1 ). The zonally symmetric tide has no phase progression with height, in contrast to the 
clear westward phase progression of the T 2am and T 2pm fields. The MarsWRF dust storm simulation does show 
an amplification of the zonally symmetric tide during the storm (increasing from 1 K to 2-3 K); however, the 
amplitude is less than half the amplitude of the simulated DW2 during most of L s = 200°-240° at low and 
middle latitudes (Figures 6a and 6d). Figure 6a shows the MarsWRF simulation of the 2001 global-scale dust 
storm mimicking this amplification with a rapid increase in DW2 amplitude at an approximate height of 
45 km from around 1 K to more than 6 K immediately following global expansion of the dust storm at 
4 = 205°, with a smaller but analogous increase in amplitude during the first phase of dust lifting near 
4 = 180°. Using a Fourier analysis of the MarsWRF temperature fields in longitude and time, the direct DW2 
tidal amplitude is seen in Figures 5b, 5c, and 6a-6c, without the aliasing implicit in the TES T-difference and T- 
average fields. The MarsWRF dust storm simulation clearly shows the strong amplification of DW2 during the 
storm (Figure 5b), exceeding 10 Kat high altitudes, compared to the very weak amplitudes seen in the control 
simulation (Figure 5c). Equatorial DW2 amplitudes in the control simulation remain less than 2 K for the entire 
simulation (Figures 5c and 6b). By 4 = 230-240° the amplitude of the observed m = 1 T-difference field 
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(o) WRF-2001 Dust Storm DW2 Amplitude (b) WRF-Control DW2 Amplitude 
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Figure 6. Simulated tidal amplitude at approximately 45 km altitude for the (a, d, g, j) WRF global dust storm simulation, 
(b, e, h, k) control simulation, and (c, f, i, I) generic dust storm simulation for four nonmigrating tides. Tides plotted are 
the (Figures 6a-6c) diurnal westward wavenumber 2 tide, (Figures 6d-6f) zonally symmetric tide, (Figures 6g-6i) diurnal 
wavenumber 1 Kelvin wave, and (Figures 6j— 61) diurnal westward wavenumber 3 tide. 

decreases by more than 50% near the equator (Figure 5d), and this pattern is duplicated in the MarsWRF dust 
storm simulation (Figures 5e and 6a). The control simulation continues to show amplitudes that are weaker 
still (Figures 5f and 6b). 

The amplification of DW2 is also seen in T 15 observations (Figure 7a), although to a lesser degree due to the 
depth weighting of the T 15 observations smoothing the amplitudes of this tide. Still, a clear enhancement and 
near 180° phase shift is seen associated with the initial dust lifting near FHellas around L s = 185°. Amplitudes 
remain higher through L s = 220° before a rapid weakening thereafter. The nearly 1 80° phase shifts indicated 
in DW2 in Figure 7a between L s = 185°-220° are represented in the MarsWRF dust storm simulation as well. 

The reason for this intense amplification of DW2 is the longitudinal structure present in the dust mixing ratios 
(and hence the thermotidal heating) during the storm. Clancy et al. [201 0] and Cuzewich et al. [201 3a] both 
describe the large wavenumber 1 pattern in dust mixing ratio seen at high altitudes (near 60 km) during the 
storm, and it is reproduced in Figures Ic-lf. The response to this intense wavenumber 1 forcing modulation is 
the strengthened DW2. To test this, we conducted a "generic" dust storm simulation without the vertical and 
horizontal structure in dust mixing ratio that TES observed with global optical depths set to 2 at a reference 
pressure of 610 Pa. Thus, this simulation has none of the zonal structure in dust mixing ratio exhibited in the 
2001 global-scale dust storm while possessing comparable total atmospheric dust loading. Indeed, with 
comparably high global dust opacity but without the wavenumber 1 pattern in tropical high-altitude dust 
mixing ratio, the DW2 tide is weak and very similar in amplitude to that seen in the control simulation (Figure 6c). 
The generic dust storm simulation suppresses most of the nonmigrating tides (Figures 6c, 6f, and 61) except for 
DK1 (Figure 6i) which tidal theory suggests should be relatively invariant with changing atmospheric dust 
loadings [e.g., Guzewich et al., 2013b] and is indeed relatively unchanged in the generic dust storm scenario. 

Not only does DW2 dominate the observed nonmigrating tide field during the L s = 200°-210° period (e.g., 
Figure 3), but DK1 seems to weaken, confirming the result of Lewis and Barker [2005] and Wilson et al. [2008a, 
2008b] in surface pressure analysis. In the MarsWRF control simulation, the mean equatorial (10°S-10°N) 
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Figure 7. Fourier fits in longitude over time in MY25 of the 10°S-10°N TES T 15 T-difference temperatures for (a) wavenumber 
1, (b) wavenumber 2, and (c) wavenumber 3. 


amplitude is 1 .5-2 K throughout the year, with the minimum occurring during the dustier near-perihelion 
portion of the year (Figure 6h). In the dust storm simulation, the amplitude is approximately 0.5 K weaker 
during and after the global dust storm (Figure 6g). This change is also reflected in the observed m = 2T- 
difference field (Figure 8) where the equatorial amplitude decreases from >4 K before the storm and in 
nondust storm years (Figures 8a and 8b) to 1-2 K in MY25 L s = 200°-210° (Figures 8d and 8e) before increasing 
modestly near 4 = 220° (Figure 8f). To confirm this weakening, we examined the TES T 15 observations. 
Figure 7b shows the m = 2 component of the T-difference temperature, which is representative of DK1 and 
confirms this slight weakening (and phase shift) after initial dust lifting at L s = 185°. DK1 remains relatively 
weak through L s = 210° when it begins to strengthen in the T 15 observations. The strengthening after 
4 = 210° seen in both types of TES observations makes sense in light of the state of the storm at that point, 
when dust lifting was occurring in both hemispheres (near Solis Planum and FHellas Basin) providing a s = 2 
component to the thermotidal forcing. DK1 is forced by the interaction of the migrating diurnal tide with the 
s = 2 component of the topography and thus an additional s = 2 component to the thermotidal forcing should 
enhance DK1 at this phase of the storm. As the dust lifting ceased after 4 = 220°, this additional forcing was 
removed and both the T 15 and limb-sounding observations show DK1 weakening yet again. 

In addition to DK1, the westward diurnal wave 3 nonmigrating tide (DW3) also contributes to the m = 2T- 
difference field, and both MarsWRF simulations produce comparable amplitudes between DW3 and DK1 
(Figures 6g, 6h, 6j, and 6k). DW3 has a much shorter vertical wavelength than DK1 . Profiles of the phase of the 
observed m = 2 T-difference field (not shown) are a means to distinguish the influence of these two tides on 
this aliased satellite-relative field. Prior to the storm, the phase of the observed m = 2 T-difference field is 
roughly constant with height near the equator, indicative of the dominance of DK1 over this field. Flowever, at 
4 = 2 1 0° there is clear westward propagation of phase with height from 40 to 60 km. By 4 = 240°, the 
observed phase looks similar to the prestorm phase diagrams, with roughly constant phase with height in the 
equatorial latitudes. We interpret this as likely brief amplification of DW3 near 4 = 200°-220° (replicated in the 
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Figure 8. Amplitude of the observed limb-sounding TES m = 2 T-difference field for (a) MY24 4 = 180°, (b) MY24 4 = 200°, 
(c) MY25 i s = 190°, (d) MY25 4 = 200°, (e) MY25 4 = 210°, and (f) MY25 4 = 220°. This represents DK1. 


MarsWRF simulation in Figure 6j) to dominate the observed m = 2 T-difference field before it weakens, and 
the two tides are roughly comparable in amplitude for the remainder of the storm until 4 = 240°. Flowever, 
the observed amplification of DW3 may simply be part of its seasonal cycle (as seen in the control MarsWRF 
simulations, Figure 6k), which was only apparent because of the weakening of DK1. Given the strong 
wavenumber 1 pattern to the thermotidal forcing (as exhibited in the pattern of high-altitude dust mixing 
ratio), we would not expect DW3 to be more strongly forced than is typical at this time of year. 

The observed m = 3 T-difference field, representing a combination of the wavenumber 2 diurnal Kelvin wave 
(DK2) and westward diurnal wavenumber 4 tide (DW4) (see Table 1), is similarly weakened during the storm. 
This is indicated in the T 15 observations in Figure 7c. Equatorial amplitudes of 4-5 K are present in the limb- 
sounding observations at 4 = 1 80°, which weaken to 1 -2 K by 4 = 200°. The amplitudes increase again after 
4 = 210° and return to 4-5 K by 4 = 250°. The MarsWRF dust storm simulation portrays a dramatic weakening 
of both the DK2 and DW4 amplitudes during the storm relative to the control simulation (by more than 50%, 
not shown). This further confirms the trend of weakening nonmigrating tides during the storm, with the 
exception of DW2 and possibly DW3. 


3.3. Stationary Waves 

Stationary planetary waves are forced by zonal wind flow over topography. On Mars, the dramatic (high- 
amplitude) topography forces wavenumber 1 and 2 (s = 1 and s = 2, respectively) stationary waves in both 
hemispheres, with the former dominating the southern hemisphere and latter the northern hemisphere due 
to the predominant topography in each [Banfield et a!., 2003]. The observed s = 1 stationary wave 
(represented in the m= 1 T-average field) is significantly amplified in the northern hemisphere during the 
storm. At 4 = 170°, very weak amplitudes are present in the m = 1 T-average field, consistent with the 
MarsWRF dust storm and control simulations. As northern hemisphere fall begins, the s= 1 stationary wave is 
expected to strengthen in the north, but the development of the global-scale dust storm enhances this 
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Figure 9. (a) MY25 L s = 210° amplitude of the TES limb-sounding observed m = 1 T-average field, (b) stationary wave 1 amplitude 
from the MarsWRF global dust storm simulation, and (c) stationary wave 1 amplitude from the MarsWRF control simulation, (d-f) 
The same fields for MY25 L s = 270°. 


process. By L s = 180°, a 6-7 K amplitude maximum is observed near 30-50 km altitude at 50°-70°N (not 
shown). This wave response continues to strengthen through L s = 210°-220°, when amplitudes reach 1 1-13 K 
(within regions of retrieved data, Figure 9a). A similar response is seen in the MarsWRF dust storm simulation, 
although in the GCM the amplitude max occurs at higher altitudes (50-70 km, Figure 9b). By contrast, the 
control MarsWRF simulation shows much weaker amplitudes (Figure 9c). Interestingly, the s= 1 stationary 
wave response is favored, while the s = 2 stationary wave remains weak. This is in contrast to the usual 
stationary wave behavior in the northern hemisphere, where the s = 2 stationary wave is favored due to the 
predominantly s = 2 topography of the northern lowlands [Banfield et al., 2003; Guzewich et al., 2012]. The 
MarsWRF dust storm simulation also produces a strong s = 2 stationary wave response in the tropics, but this 
is not observed in the TES data. 

The s= 1 stationary wave in the northern hemisphere may be a response to a westward propagating 
equatorial wavenumber 1 structure that appears during the 2001 global-scale dust storm as indicated by TES 
limb-sounding observations of dust mixing ratio (Figures Ic-lf) [ Clancy et al., 2010; Guzewich et al., 2013a]. 
Interestingly, the observed phase of the s= 1 stationary wave near the amplitude maximum in the northern 
hemisphere is not constant during the storm (as would be expected without the storm) [e.g., Guzewich et al., 
2012] but follows the high-altitude dust mixing ratio pattern with time before "normalizing" to the expected 
value of ~30-50°E after L s = 270°. Due to this steady change in phase, we will refer to this wave as "quasi- 
stationary." This may in fact also be an indication of a slow-period traveling wave, similar to that described by 
Wilson et al. [2002] and Banfield etal. [2004], although eastward traveling waves would normally be expected 
in this season. This change in phase is demonstrated in Figure 1 0, where the 60 km altitude TES dust mixing 
ratios from Figure 1 are plotted for L s = 190-210°. Figure 10 (solid black line) indicates the phase of the 
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Figure 10. TES 2 P.M. dust mixing ratio (Ar/mb) at 60 km altitude for MY25 (a) L s = 190°, (b) L s = 200°, and (c) 4 = 210°. Solid black 
line indicates the phase (in east longitude) of the maximum amplitude of the northern hemisphere wavenumber 1 stationary wave. 
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maximum amplitude of the northern hemisphere s = 1 quasi-stationary wave at each time step. During the 
L s = 190°-220° period, when the equatorial wavenumber 1 structure is most obvious in the dust mixing ratio 
plots at high altitude (50-60 km), the phase of the northern hemisphere s = 1 quasi-stationary wave (represented 
by the m = 1 T-average field) directly corresponds to the leading (western) edge of the higher dust mixing ratios 
and the eastern edge of the "dust-free corridor" described by Clancy et al. [201 0]. Wilson et al. [2002] suggested 
that tropical inertial instability maintains an eastward propagating slow s = 1 traveling wave following the 
results of Salby and Callaghan [2001]. Although we cannot directly test for observed tropical inertial instability in 
the TES data, we can do so in the MarsWRF simulations. We find a large region of inertial instability (as 
manifested by negative absolute vorticity in the northern hemisphere) is present in the low latitudes (0-1 0°N) 
of the dust storm simulation during the L s = 1 90°-220° period that is absent in the control simulations (not 
shown). This suggests that during the period of highest dust opacity, tropical inertial instabilities did at least 
maintain the quasi-stationary or slow-traveling wave seen in the northern hemisphere during the storm. 

Even after the storm decayed at i 5 = 270°, the observed TES s= 1 quasi-stationary wave amplitude remains 
relatively weak (Figure 9d) which is generally reflected in the MarsWRF dust storm simulation (Figure 9e) 
relative to a "normal" year seen in the control simulation (Figure 9f) with much larger wave amplitude. 

4. Conclusions 

The 2001 (MY25) global-scale dust storm initiated dramatic changes to the wave structure of the Martian 
atmosphere. Using temperature retrievals and brightness temperatures from TES observations and supporting 
information from a MarsWRF general circulation model simulation of the dust storm, we have identified several 
significant and unique changes to the thermal tides and stationary waves forced by the dust storm. 

As expected from previous studies [e.g., Wilson and Hamilton, 1996], the vertically trapped extratropical 
modes (in the summer hemisphere) of the diurnal migrating tide are strongly forced by the high-column dust 
opacity during the storm, with amplitudes exceeding 20 K in the southern hemisphere and 10 K in the 
northern hemisphere. The equatorially confined gravest Flough mode of the diurnal migrating tide is 
suppressed to undetectable levels and does not reemerge until after Z. s = 270°. 

Nearly all nonmigrating tides, including the diurnal Kelvin waves (DK1 and DK2), are repressed during the 
early portion of the storm's global extent at F s = 200°-210°. The notable exception to this pattern being the 
diurnal westward propagating wavenumber 2 tide and possibly the diurnal westward propagating 
wavenumber 3 tide. This highlights the importance of the zonal distribution of aerosol forcing. Theoretical 
predictions and some previous work have shown that the wavenumber 1 diurnal Kelvin wave's temperature 
amplitude should remain relatively invariant to zonally uniform changing atmospheric dust loading [e.g., 
Guzewich et al., 201 3b]. The changing pattern of dust lifting centers during the storm and the predominantly 
wavenumber 1 forcing by airborne dust invalidated that zonally uniform criterion. As a result of this 
wavenumber 1 dust forcing, the diurnal westward propagating wavenumber 2 tide, which appears in the 
m = 1 T-difference field in the satellite-relative perspective of TES, was strengthened. The amplitude of this 
tide increases to 4-8 K, with a similar amplification present in the MarsWRF dust storm simulation. 

During the dust storm, dust at high altitude appeared to be tied to a slowly westward propagating wavenumber 
1 wave [Clancy et al., 2010; Guzewich et al., 2013a]. This appears to have forced and/or maintained a strong 
(> 1 0 K) wavenumber 1 northern hemisphere quasi-stationary (or slow-traveling) wave. The phase of this wave is 
consistent with its maximum located at the leading (western) edge of the high dust mixing ratios. At the same 
time, the wavenumber 2 stationary wave is rather weak in the northern hemisphere, contrary to typical 
conditions during this season of the year [Guzewich et al., 2012]. 

This study demonstrates the radical alteration by a global-scale dust storm to the wave forcing mechanisms 
normally present on Mars and provides a broader planet-wide view of the atmospheric dynamics during the 
storm. This adds confirmation to the coupling of waves between the tropics and higher latitudes (as described 
in Wilson et al. [2002]) through the link between the equatorial wave apparent in dust mixing ratio 
observations and associated inertial instability with the unusually strong high-latitude quasi-stationary wave. 
The demonstrated importance of the zonal structure in the dust mixing ratios confirms suggestions from 
previous studies that wave forcing in the atmosphere is more than a simplistic function of total atmospheric 
optical depth, which might be a more reasonable assumption if only the migrating tides were examined. 
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